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ABSTRACT 

We report the results of a low frequency radio variability and slow transient search using archival 
observations from the Very Long Array. We selected six 325 MHz radio observations from the spring of 
2006, each centered on the Spitzer-Space- Telescope Wide-area Infrared Extragalactic Survey (SWIRE) 
Deep Field: 1046+59. Observations were spaced between 1 day to 3 months, with a typical single- 
epoch peak flux sensitivity below 0.2 mJy beam _1 near the field pointing center. We describe the 
observation parameters, data post-processing, and search methodology used to identify variable and 
transient emission. Our search revealed multiple variable sources and the presence of one, day-scale 
transient event with no apparent astronomical counterpart. This detection implies a transient rate of 
1+1 event per 6.5 deg 2 per 72 observing hours in the direction of 1046+59 and an isotropic transient 
surface density £ = 0.12 deg~ 2 at 95% confidence for sources with average peak flux density higher 
than 2.1 mJy over 12 hr. 

Subject headings: methods: observational — radio continuum: general — stars: oscillations — stars: 
variables: other 



1. INTRODUCTION 

A long list of sources may prese nt variable and tran- 
sient behavior, from stellar flares dBastian et al. | | 1998l: 



varia ble and transien t radio s ources (iHvman et all 
20051: lLorimer et all 120071: IBower et al I 12007 ' 



Osten et alll2005h, spinning neutron stars (jHewish et all 
19681: ICamilo et all \mi&i. and gamma ra y bursts 
(jKlebesadel et all 119731 : iDessenne et all 11996ft . to yet- 
detected extra-solar plane ts dZarkal 119981 ) and gravita- 
tiona l wave counterparts (jBlancheti 2002; A bbott et all 
2009). Radio observations serve an important role in the 
investigation of these objects and their variable behavior. 
Measurements at radio wavelengths uniquely probe mag- 
netic field topology and non-thermal processes. This pro- 
vides key information for understanding the acceleration 
mechanisms responsible for intense, time-variable emis- 
sion. Past exploration of radio variability has been signif- 
icantly limited in frequency range, time resolution, and 
spatial coverage. Recent improvements in low frequency 
radio spectrometer sensitivity, radio frequency interfer- 
ence (RFI) mitigation, and wide-field imaging techniques 
now place radio observations at the fore-front for the de- 
tection and subsequent understanding of these dynamic 
objects. 

Early blind searches have already uncovered multiple 
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Burke-Spolaor fc Bailed 
However, a majority of the observations were performed 
at gigahertz (cm to mm wavelength) frequencies, leaving 
the sub-gigahertz (meter wavelength) sky vastly unex- 
plored. We present the results from a sub-gigahertz 
variability and transient search using data from the 
Very Large Array (VLA) archive. We describe the 
observation and analysis of 6x12 hr observations at 
325 MHz (0.9 m, 6.5 deg 2 field of view) with a common 
pointing center and temporal spacing between 1 day 
and 3 months. We report trends in the variability of 
approximately 950 field sources, highlighting multiple 
sources with significant flux variability. Further, we 
report the discovery of one transient event unassociated 
with a previously cataloged source and discuss the 
implications of our detection on existing transient limits. 

2. OBSERVATIONS AND ANALYSIS 

2.1. Radio Observations 

We analyzed six 325 MHz radio observations from 
the VLA archive centered on a = 10 ,I 46 m , S = +59 d 
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TABLE 1 

Summary of radio observations listing the 2006 start date, observation time range (UTC), bandpass calibrator, amplitude 

AND PHASE CALIBRATOR, NUMBER OF UV VISIBILITIES AFTER FLAGGING, AND RMS RADIO FLUX (mJy BEAM -1 ) AT THE POINTING CENTER. 



Start Date 


Time Range 


BP 


A&P Cal 


Visibilities 


RMS 


02/19 


02:15:59 - 14:13:39 


0137+331 


0542+498 


1660567 


0.179 


02/27 


01:45:20 - 13:42:19 


0137+331 


0137+331 


2076319 


0.171 


03/03 


01:30:29 - 03:26:29 


0137+331 


0542+498 


1622207 


0.184 


03/04 


01:25:49 - 13:22:40 


0137+331 


0542+498 


1379554 


0.190 


03/10 


01:02:59 - 12:58:59 


0137+331 


0542+498 


1725216 


0.174 


05/16 


20:36:09 - 08:31:34 


0137+331 


0542+498 


1395740 


0.192 



recorded between February 19th, 2006 and May 17th, 
2006 (see Tab. [T]). The data was originally part 
of a larger effort by F. Owen to study the multi- 
wavele ngth properties of the general radio source popu- 
lation (lOwen et al.ll2009i) . Individual field measurements 
have nearly identical observation parameters, providing 
an ideal dataset to search for source variability and tran- 
sient emission. 

Each observation lasted 12 hours and utilized two 
6.25 MHz bands centered at 327.5 MHz and 321.6 MHz. 
Approximately 11 hours were allocated for measurements 
of the target field, with the remaining time allocated for 
calibrator observations. Data was collected in pseudo- 
continuum mode, with 15 spectral channels requested per 
band to mitigate radio frequency interference (RFI) and 
reduce bandwidth smearing of bright out-of-field confus- 
ing sources. The radio array consisted of 23, 25 m anten- 
nas in A- Configuration during each observatiorQ (max- 
imum baseline ~ 35 km), resulting in an approximate 
6" x 5" angular resolution and ~ 6.5 deg field of view 
(based on a 2.9 deg full- width-half-power primary beam, 
see Sec. US]). 

2.2. Field Calibration and Imaging 

To ensure uniform processing of the radio data, we de- 
veloped an automated reduction pipeline based on stan- 
dard routines from the AIPS and 0&it0 software pack- 
ages. Each observation was processed as follows: 

1. Initial data flagging was performed to eliminate two 
channels with known spectral RFI and all UV sam- 
ples stronger than 50 Jy. 

2. UV data was calibrated using standard AIPS tasks. 
The receiver bandpass correction was determined 
using calibration source 0137+331. The radio flux 
scale and antenna phase offsets were determined 
using a single 5 minute scan of 0542+498 (assumed 
46.0 Jy at 327.5 MHz, 46.2 Jy at 321.6 MHz) for all 
epochs except February 27. For this epoch, all cal- 
ibration was performed using 0137+331 (assumed 
41.5 Jy,41.8 Jy) due to substantial RFI present dur- 
ing the 0542+498 scan. 

1 Array re-configuration began before the final observation, how- 
ever, only the sub-set of antennas and locations common to all 
epochs were chosen for our analysis. 

2 Astronomical Image Processing System, release 31DEC10 

3 Obit is developed and maintained by Bill Cotton at The Na- 
tional Radio Astronomy Observatory in Charlottesville, Virginia, 
USA and is made available under the GNU General Public License, 
version 1.1.269-6-64b. 



3. An additional antenna phase correction was per- 
formed by referencing a field source model based 
on the Westerbork Northern Sky Survey (WENSS, 
IRengelink et alj|l997ft . 

4. Each field was imaged and CLEANed to a mini- 
mum flux level of 1 mJy. Field sources were sub- 
tracted to create a residual UV dataset. A second 
(and more extensive) round of RFI excision was 
performed by comparing the original UV data to 
the residual, using Obit tasks LowFRFI and Aut- 
oFlag. 

5. The fully flagged UV data was re-imaged using 
Obit task Imager which is similar to AIPS IMAGR 
but can also execute self-calibration loops and 
auto-identify bright confusing sources. We pre- 
formed two rounds of phase-only self-calibration 
and an addition round of amplitude and phase self- 
calibration for each observation. Final RMS am- 
plitude and phase variations were generally below 
3% and 0.5 deg, respectively. Field images were 
made with a 3 deg diameter field-of-view (slightly 
larger than the primary beam FWHP) and 1" pix- 
els. Radio sources were automatically windowed 
and CLEANed to a level of 0.1 mJy. 

2.3. Source Identification and Filtering 

Sources whose peak flux density exceeded 3 times the 
RMS flux density were cataloged using Obit task Fnd- 
Sou. Individual catalogs were created for each epoch. 
Source characteristics were determined using a two- 
dimensional gaussian fit, providing an estimate of the 
center position, shape, peak and integrated flux density, 
local background RMS (within 300"), and amount of 
background slope or curvature. The catalogs were au- 
tomatically filtered to exclude objects smaller than the 
synthesized beam. Highly extended sources (i.e. diame- 
ter > 30" ) were not well modeled by a single 2D gaussian 
and were also excluded. All flux estimates were corrected 
for primary beam attenuation. 

Sources from the six individual catalogs were matched 
across epoch by center position. A 5" matching radii 
(approximately one synthesized beam width) was used 
to account for fit error and/or any potential ionospheric 
refraction. Sources positioned within 5" were further dis- 
criminated using the peak flux density. The measured 
variation in source position (Aa, AS) between the mean 

4 Global 3<r before primary beam correction. We anticipate a 
large number of " accidental" detections at this threshold. A second 
threshold is applied later, based on Gaussian noise expectations. 
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position and the position observed for each epoch is plot- 
ted in Figure [TJ Changes in source position were gen- 
erally less than 1" for all locations inside the primary 
beam. All displacements greater than 2" corresponded 
to inadequate automated fitting of extended sources with 
complex structure (see Sec. 13.11 for more discussion on 
the behavior of resolved and unresolved sources). Fig- 
ure [2] displays the average single epoch RMS (la) radio 
flux versus distance from the pointing center for sources 
persistent in all observations. The RMS measurements 
ranged from 0.19 mJy beam -1 to 0.41 mJy beam -1 and 
roughly trace primary beam attenuation profile. RMS 
values indicate a ~ 1.45 deg half-power radius and no 
discernible beam asymmetry. 

The combined catalog was divided into portions rep- 
resenting (Group 1) sources detected in all 6 epochs, 
(Group 2) sources detected intermittently, i.e. sources 
found in multiple but not all epochs, and (Group 3) 
sources detected in only one epoch. To help identify de- 
tections in Group 2 and Group 3 which are associated 
with low signal-to-noise field sources (sources in a given 
epoch with peak flux below the initial global 3cr thresh- 
old), we created and cataloged a composite image made 
from the 6 individual observations. This image had a 
RMS < 0.1 mJy beam -1 near the pointing center and 
contained over 1300 sources with integrated flux density 
exceeding 0.5 mJy. The number of sources and the inte- 
grated flux density values in our composite catalog were 
consistent with nu mber and flux mea surements of the 
same field made bv lOwen et al.1 (|2009f ). Any single and 
intermittent detections coincident (5" search radii) with 
sources in the composite catalog were reevaluated to de- 
termine their properties for all 6 epochs. Approximately 
90% of the intermittent detections were reclassified and 
added to the list of persistent field sources (i.e. Group 1) 
during this process. No single epoch detections were re- 
classified during this process. All sources in Group 1 were 
searched for evidence of variability (see Sec. I2.4[) . The 
remaining source detections from Group 2 and Group 3 
were filtered by signal-to-noise ratio to identify statisti- 
cally significant transient events. This filter is described 
in Section [231 

2.4. Measuring Source Variability 

We searched the list of field sources identified in all ob- 
serving epochs (Group 1) for variability in the measured 
total flux density. We quantify the flux variation using 
two methods, calculating the (1) modulation index m, 
i.e. standard deviation divid ed by the mean as used by 
iGaensler fc Hunsteadl (|2000ft . and (2) a linear fit span- 
ning the 6 epochs. Moreover, we checked all sources for 
any variation in shape. Changes in source shape (A0) 
were quantified by first computing the standard devia- 
tion in the source major and minor axis from their mean 
value, normalizing each result by the mean value, and 
finally summing in quadrature. 

We assume that the flux density for a majority of 
the radio sources is stable over timescales longer than 
3 months. However, source measurements from the indi- 
vidual epochs appeared uniformly offset from their mean 
value, presumably caused by long-duration amplitude 
changes in the primary amplitude calibrator. To account 
for this effect, we normalized each observation to the 
composite image catalog described in Section [231 The 
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Fig. 1. — Variation in the measured position (Aa, AS) of field 
sources identified in all observing epochs. Compact sources were 
generally located within 1" with larger separations typically asso- 
ciated with inconsistent fitting of AGN. 
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Fig. 2. — Average RMS radio flux versus radial distance from 
the field center near sources detected in all epochs. Values are 
determined via histograms of the radio flux measurements within 
300" of each source. The RMS measurements trace the attenuation 
of the primary beam, indicating a half-power radius of ~ 1.45 deg. 

gain corrections (one for each epoch) were determined by 
minimizing the standard deviation in total flux of the 100 
brightest field sources from their composite value. The 
resulting gain corrections ranged between 2% (February 
27) and 14% (May 16). These corrections are consistent 
with VLA reported flux density variations for calibrator 
sources 0137+331 and 0542+498 during our observation 
period. Further, gain corrections did not appear spatially 
correlated. Performing the above normalization process 
for sub-regions inside the field of view produced similar 
results. 

2.5. Transient Source Detection 

We searched the list of intermittent (Group 2) and 
single-epoch (Group 3) detections for previously unde- 
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Fig. 3. — (full) 325 MHz radio image centered on a = 10' i 46 m , 6 = +59 d observed February 19, 2006. The image has been corrected for 
primary beam attenuation, (inset) Flux density histogram. The image noise distribution is approximately Gaussian. 



tected transient events. A vast majority of the "candi- 
date" transients were due to noise fluctuations or sys- 
tematic effects corresponding to bright field sources. As- 
suming the image noise can be modeled using Gaussian 
statistics (cf. Fig. [3] inset for verification), the expected 
number of transient detections N above a given threshold 
n a in a given field can be estimated as 

N = n„- erfc{n a /V2), (1) 

where erfc is the complementary error function and n s is 
the number of independent samples 

Tlobs X Q Q bs 



given by the number of epochs n bs times the observable 
area for each epoch £l b s divided by the area of the syn- 
thesized beam flf,. Note that, in general, the number of 
independent samples n s is also dependent on the number 



of UV visibilities. For example, an image made from only 
the longest observation baseline will have no indepen- 
dent samples. However, as we are primarily focused on 
long duration transients and our VLA observations gen- 
erate a large number of uniformly spaced visibilities per 
hour, we have omitted this dependency from our analysis. 
Such consideration may be required to analyze smaller 
timescales or when using other radio arrays. Further, the 
maximum observable area varies slightly with threshold 
choice due to diminishing beam sensitivity away from the 
pointing center (cf. Fig. [2]). In practice, one chooses a 
flux threshold such that the predicted number of " acci- 
dental" detections is much less than one. For this obser- 
vation, n s is approximately 2 x 10 7 using a search area 
which extends to the beam half-power points and N is 
unity for a signal-to-noise threshold n„ ~ 5.5. We de- 
fine the signal-to-noise as the ratio of the source fitted 
peak to the local RMS noise. Flux measurements larger 
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than 5.5a are then either due to real sources or indicate 
systematic effects. 

Accidental detections associated with systematic ef- 
fects were mitigated by excluding candidates within 10' 
of source with median peak flux exceeding 100 mJy 
beam -1 . A grey-scale image of the 1046+59 field is 
shown in Figure [31 The image has been corrected to re- 
flect the primary beam attenuation profile. Circles mark 
the 12 sources brighter than 100 mJy beam~ 1 and the 
corresponding 10' regions used to filter accidental detec- 
tions caused by strong residual side-lobes. The excluded 
regions account for an approximate 5% reduction in the 
transient search area fl b s - 

Sensitivity variation across the primary antenna beam 
(cf. Fig. [2]) manifests as a trade-off between the ob- 
serving area ft b s and the desired detection threshold 
n a . Figure [4] illustrates the measured number of acciden- 
tal detections versus distance from field pointing center. 
Curved lines indicate the expected number of acciden- 
tal detections N for different thresholds n a> assuming 
Gaussian field noise (cf. Fig. [3]inset for verification) and 
the measured primary beam attenuation profile. For this 
search, threshold of 5.5tr will detect significant events at 
any location inside the primary beam and maximize il bs- 
Detections exceeding 4.5<r are significant at all locations 
within 0.1 deg (~ 60-70 synthesized beams) of the point- 
ing center. However, while the decreased threshold pro- 
vides a factor of 2.5 better flux density sensitivity, the 
observing area for 4.5c events decreases by more than a 
factor of 200. 

To identify statistically significant transient sources, 
we filtered the single-epoch detections to exclude sources 
with peak flux density below a 5.5cr, maximizing fi f, s . 
We filtered the intermittent detection list at a slightly 
lower threshold (5a, equal to the minimum flux threshold 
for a single epoch) to provide more sensitivity to tran- 
sients which occur between observations. To further eval- 
uate the candidates, we required any positive transient 
event to be centered on a 300" diameter region free of 
spurious noise fluctuations at a 90% expectation level, 
i.e. a background noise threshold of 4.1er (n s ~ 2500, 
N ~ 0.1 for n a — 4.1). This is an important final step. 
Regions failing this criteria display non- Gaussian noise 
properties (likely due to systematic imaging effects) and 
the significance of any transient can not be determined 
by a single threshold test. Sources which remained after 
filtering were marked for visual inspection. 

2.6. Transient Detection Rate 

The number of anticipated transient events N t is the 
product of the total observing area n a b s x f2 Q & s times the 
isotropic tra nsient surface density E. Measurements at 
843 MHz bv lBannister et all (|2011f ) suggest that £ may 
be as high as 1.3 x 10~ 2 dcg~ 2 for long duration (day 
to year) transients stronger than 8 mJy. Assuming flat 
spectrum sources, we expect ~ 1 transient detection from 
this search. However, there is still little known about 
the transient surface density frequency dependence and 
if scaling £ in such a manner is valid. 

The probability of detecting N transients from an ex- 
pected number of events Nt is given as 




1.0 0.5 1.0 1.5 

Distance from Pointing Center (deg) 

Fig. 4. — Number of accidental detections versus distance from 
the pointing center. Histograms indicate measured values. Solid 
curves indicate detection estimates for different flux thresholds n CT , 
assuming a Gaussian field noise distribution and incorporating the 
primary beam attenuation profile. A single transient detection is 
statistically significant only when the predicted number of 'acci- 
dental' detections at the same flux or SNR threshold is less than 
unity (dashed line). 

assuming a Poisson distribution of sources. If Nt is un- 
known (i.e. an unknown surface density), the measured 
number of transient detections can be used to establish 
an experimental upper limit for a set confidence level. In 
the limiting case of no detections (N = 0), the expected 
number of events Nt is approximately 3 at 95% condi- 
dence (P — 0.05) and the upper limit to the isotropic 
transient surface density can be expressed as 



£ < 



3.0 



n b s x n o b s (n a ) 



(3) 



P(N-Nt) 



Np—Nt 



Nfe 



(2) 



Qobs(n<j) indicates that, in general, the observable area 
depends on the flux density threshold. 

3. SEARCH RESULTS 

3.1. Field Source Variability 

Approximately 950 field sources were persistent in all 
six observations above an integrated flux density of ~ 
0.6 mJy (3a at the field center). Each of theses sources 
was searched for evidence of variability in the total flux 
density using the methods outlined in Section 12.41 We 
find that the random variability in field sources sources 
is low and changes in a relatively uniform way. Most de- 
tections have a modulation index on the order of 10% for 
SNR > 20 and below 5% for SNR > 50. The observed 
source variation is directly comparable to the fractional 
variations in the Gaussian major and minor axis fit pa- 
rameters. 

We observed excess variability in a small number of 
high SNR detections (~ 25 total sources, SNR > 50). 
For a majority of these detection (> 80%), the mea- 
sured modulation index appears associated with inade- 
quate automated fitting of resolved sources with complex 
structure. In these specific cases, fit parameters (center 
position, major/minor axis, etc.) varied by as much as 
50% from the mean value and caused anomalous posi- 
tion and flux measurements. In the remaining cases, ex- 
cess variability appears consistent with interstellar scin- 
tillation of a compact components in extragalactic radio 
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sources. Meter- wavelength observations of ~ 10 mas ra- 
dio sources often display amplitude fluctuations ranging 
from 3% t o 10% over timescales of days to mo nths, re- 
spectively (|Blandford et al.lll986t lRickettHl986T >. 

Table [2] lists sources which have a modulation index 
m greater than that observed for field sources of similar 
SNR, but are not easily explained by variations in source 
sha pe A8 caused by fit irregularities (m > A8, see Sec. 
[2T4|) . Redshift measurements from the NASA/IPAC Ex- 
tragalactic Database (NED) are also given, if available. 
Figure [5] displays the 3-month flux variation (standard 
deviation divided by the mean) versus the signal-to-noise 
ratio (fitted peak divided by RMS) for all persistent field 
sources. Circles indicate measurements corresponding to 
unresolved sources, while squares indicate measurements 
corresponding to resolved sources. The brightest sources 
in the field are predominantly radio galaxies and are of- 
ten resolved with meter- wavelength VLA observations. 
Sources with obvious fit irregularities have been excluded 
from Figure [5j roughly characterized by shape variations 
greater than 15% for sources with SNR > 50. 

In addition to measurements of random variability, 
we also monitored each persistent field source for lin- 
ear trends, i.e. sources with significant brightening or 
dimming over the 3 month observation period. The two 
sources with the largest flux density trends are high- 
lighted below. 

J 1047 19.1+ 5821 1 7 

J104719.1+582117 is optically identified as a 
quasar with an appr oximate distance of 350 Mpc 
(|Abazaiian et a l. 2009). We observed a steady decrease 
in the flux density of J104719.1+582117 across the 6 
observing epochs (see Fig. [51 top). The measurements 
ranged from 47.7 ± 0.39 mJy on February 19th to 
35.1 ± 0.39 mJy on May 17th, a 26% decline in flux and 
a linear best-fit slope of -0.143 mJy/day. 

J104539.1+580711 

J104539. 1+580 711 is a quasar locat ed at approxi- 
mately 330 Mpc (jAbazaiian et al.l 120 09). The peak ra- 
dio flux density increased steadily over the observation 
epochs, growing from 27.1 ± 0.42 mJy on February 19th 
to 35.8 ± 0.45 mJy on May 17th. The best fit slope is 
0.094 mJy/day (cf. Fig. bottom). 
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Fig. 5. — Random variation in the 325 MHz radio flux density 
over 3 months vs. SNR (cr) for sources in the 1046+59 field. Unre- 
solved sources are plotted with open circles, while resolved sources 
are plotted with open squares. 
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Fig. 6.— Radio light curves for J104719. 1+582117 (top) and 
J104539.1+580711 (bottom) from February 19, 2006 to May 17, 
2006. Open circles indicate the integrated flux for the variable 
sources measured on each epoch. For reference, open diamonds 
indicate flux values for a nearby source of similar intensity. Dashed 
lines represent the best linear fit to the data. Uncertainty estimates 
are roughly equal to the symbol size. 



3.2. Transient Sources 

We inspected a list of 25 intermittent and 64 single- 
epoch transient candidates with peak flux density ex- 
ceeding detection thresholds of 5a, 5.5a, respectively. All 
25 intermittent detections and approximately 20 of the 
single-epoch detections were identified as weak persis- 
tent sources not recorded in our median catalog. Only 
one single-epoch source satisfied the final selection cri- 
teria given in Section 12. 5[ i.e. centered on 300" region 
free of noise fluctuations exceeding the 90% expectation 
level. The remaining candidates were "false detections" 
which were marginally detected (44 of 46 below 5.7c) 
and in 300" regions plagued with imaging artifacts from 
a nearby field source. For completeness, we have ex- 
cluded these regions from the final observation area (~ 



1% reduction to £l bs), consistent with our mitigation of 
systematic effects around the 12 brightest field sources 
(cf. Sec. 12. 3[) . More details on the one remaining source 
are presented below. 

J103916. 2+585124 

A transient radio source located at a = 10 ,I 39 m 16.2 s , 
S = +58 d 51 m 24 s was detected in the March 3rd 
epoch with peak flux density equaling 1.70 ± 0.25 mJy 
beam" 1 (6.7cr, 12 hr integration) and is centered on 300" 
region free of noise fluctuations in excess of 2.9cr. The 
source is unresolved and located 0.76 deg from the point- 
ing center (15% primary beam attenuation correction). 
J103916. 2+585124 is positioned appro ximately 5" from 
an op tical galaxy of undetermined type (|Abazaiian et al.l 
2009), but does not appear associated given the accuracy 
of the optical and radio position measurements (cf. Fig. 
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TABLE 2 

Sample variable sources in the 1046+59 field, listing the position (a, 8), Average RMS in mJy beam -1 , Average SNR, 

MODULATION INDEX m IN PERCENT, PERCENT VARIATION SOURCE SHAPE AO, AND PUBLISHED REDSHIFT Z IF KNOWN. 



Criteria 


a 


.5 


RMS 


SNR 


m 


AO 


z 




10 47 19.1 


58 21 17 


0.20 


213.30 


9.95 


5.30 


1.22 


m > 7.5 


10 45 28.3 


59 13 27 


0.18 


188.90 


10.54 


5.63 


2.31 


SNR > 50 


10 45 39.1 


58 07 11 


0.22 


115.51 


9.40 


5.27 


1.15 




10 44 50.5 


59 19 27 


0.18 


98.63 


7.69 


5.84 


0.84 




10 40 34.5 


59 54 45 


0.26 


48.31 


15.05 


3.89 




m > 10 


10 48 24.0 


58 30 27 


0.22 


38.46 


10.15 


7.24 




SNR > 20 


10 50 21.7 


58 42 57 


0.24 


27.68 


14.72 


6.97 






10 44 56.3 


59 38 02 


0.19 


21.26 


14.74 


6.07 





[7]). There are no published IR, X-ray, or GRB counter- 
parts. 

We reprocessed the March 3rd epoch in multiple ways 
to detect potential processing artifacts. We imaged the 
UV data isolating each frequency sub-band, each po- 
larization, and without self-calibration. We also im- 
aged the UV data with various pixel and facet sizes. 
The transient source was detected in all configura- 
tions. J103916. 2+585124 appears unpolarized within the 
recorded error, measuring 1.76+0.32 mJy beam -1 in right 
circular polarization (RR), 1.63+0.34 mJy beam -1 in left 
circular polarization(LL), and is undetected at a 2.5<r 
level of 0.68 mJy beam -1 in Stokes V . The source has 
no discernible frequency dispersion, with nearly equal 
intensity in each frequency sub-band (1.67 ± 0.31 mJy 
beam^at 327.5 MHz, 1.72 ± 0.31 mJy beam^at 
321.6 MHz). We did not detect J103916.2+585124 in 
images with frequency bandwidth less than 6.25 MHz. 

A plot of the peak flux density versus time is shown 
in Figure M for J103916. 2+585124 and a nearby field 
source J103946. 5+585405. Markers indicate the peak 
value from a 2D Gaussian fit at 12 hr (diamonds), 6 hr 
(squares), and 1 hr (circles) intervals. We did not de- 
tect the source on timescales shorter than 1 hr, search- 
ing for emission on temporal spacings as fine as 10 s. 
Variability is evident on sub-day timescales, with the 
source appearing weak first 6 hours, then rising by ap- 
proximately 1 mJy beam" 1 to a maximum intensity of 
2.1 + 0.29 mJy beam" 1 . The recorded signal-to-noise in- 
creases from 6.7 to 7.3 when only imaging the second half. 
J103916. 2+585124 is undetected 4 days earlier on Febru- 
ary 27th at a 2.5<r level 0.60 mJy beam -1 and 12 hours 
later on March 4th at a 2.5a level 0.63 mJy beam" 1 . 

J103738. 4+592001 AND J104304. 7+591726 

We note the presence of two events detected in the May 
16th, 2006 epoch which have a peak flux exceeding 5.5<r 
but appear in regions with non-Gaussian noise character- 
istics. J103738. 4+592001 was observed with a peak flux 
of 1.74 + 0.29 mJy beam- x (6.1cr) and J104304.7+591726 
was observed with a peak flux of 1.21 ± 0.20 mJy 
beam- 1 (6.0ff). Similar to J103916.2+585124, each 
source is visually identified in 12 hr integration images 
which isolate each observing sub-band, Stokes R and 
Stokes L, and processing without self-calibration. Nei- 
ther source has a contemporaneous optical, IR, X-Ray, 
or GRB counterpart. However, multiple noise fluctua- 
tions of comparable flux density (> 5<r, < 4.1c expected 
for 90% false detection) are observed in close proximity 



Mar 3 rd , 2006 




10 h 39 m 30 s 39 m 00 s 

IPOL 324.5 MHz 

Fig. 7.— Transient event J103916.2+585124 on March 3rd, 2006. 
The radio image scale is linear, ranging from -2 mJy bcam -1 to 
2 mJy beam -1 . The image RMS is 0.24 mJy beam -1 . White 
dotted lines outline the position of near-by galaxies m easured by 
the Sloan Digital Sky Survey (SDSS, Abazajian et aT] |2009Tl . 

to each source and the region noise could not be improved 
through reprocessing. Further, each source is undetected 
when dividing the observation into 2 x 6 hr intervals. 

4. DISCUSSION 

We detect one transient source exceeding a 5.5ct min- 
imum flux threshold (6.7tr observed) located in a region 
free of spurious noise fluctuations. The measured tran- 
sient rate in the direction of 1046+59 is then equal to 
1+1 event per 6.5 deg 2 per 72 observing hours for sources 
with average peak flux density higher than 2.1 mJy over 
12 hr. We use Equation [2] to report a corresponding 
isotropic surface density of E = 0.12 deg~ 2 95% confi- 
dence (N = l,Nt = 4.5). Figure [9] displays the isotropic 
transient surface density inferred by this search along 
with results from past observations. Multiple transient 
searches have been published with time resolution rang- 
ing from sub-millisecond to multiyear. For a basis of 
comparison, we only depict searches which are sensitive 
to day to month temporal events. All published lim- 
its have been scaled to represent the Surface Density at 
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Fig. 8. — Radio light curves for transient source 
J103916. 2+585124 and a nearby field source on March 3, 
2006. Open markers indicate the peak flux value at 12 hr 
(diamonds), 6 hr (squares), and 1 hr (circles) timescales with 
1<t error bars. Visual detections with peak flux below 2.5er are 
indicated with crosses and downward arrows. The remaining 
non-detections are drawn with a 2.5<r upper limit. 
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Equation [5] considering a minimum radio detection level S m in = 
5.5 X 0.38 mjy over 12 hours. Diagonal lines represent values for 
sources of various radii. A horizontal dashed line at 10 12 K denotes 
the approximate temperature boundary between coherent and in- 
coherent emission processes. 
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Fig. 9. — Surface density versus radio flux density for tran- 
sient searc hes with da y -scale time res olution. Lin e s indi c ate re - 
sults from IFrail et alj i2003 ) (F03), IBower et all 120071, \20W) ; 
IBower fc Saull d2011]) (B 7, B1 0, BID, IMatsumura et al] 120091) 
(M09), IBannister et at] (20TS) (Ball), and IBell et aTT MiTS) 
(Bell). Arrows indicate 95% confidence upper limits for searches 
resulting in no detections (N = 0). Estimates based on this work 
(J*) are shown in red. 



95% confidence. The IBower et ail (|2007j ) results have 
also been adjust ed to reflect the recent reanalysis by 
IFrail et all ()2011[ ). For completeness, we note the non- 
detection limit from lLazio et all ()2010l ) of S < 10~ 7 with 
flux sensitivity > 2.5 kJy at 74 MHz. This result has been 
omitted from Figure [51 

Our detection indicates a transien t surface dens i ty tha t 
is more strict than the limit from IBower et al.l (|2007| ) , 
even considering the factor of 2 decrease suggested by 
IFrail et all (|201lD . However, the frequency range in Fig- 
ure [S] (0.3 GHz to to 8.4 GHz) is greatly suppressed and 
this search was performed at frequency nearly and an or- 
der of magnitude lower lower than previous results with 
similar flux sensitivity. In general, gigahertz radio obser- 
vations typically probe non-thermal synchrotron sources 
while sub-gigahertz observations are sensitive to plasma 
effects. Different searches may be sensitive to different 



source populations. 

Transient J103916.2+585124 has no X-Ray or GRB 
counterpart, and is located approximately 5" from op- 
tical galaxy SDSS J103915.88+585127.7. The tran- 
sient emission my be explained by interstellar scin- 
tillation of a compact radio source associated with 
SDSS J103915.88+585127.7 which is below the detection 
threshold. However, the source is positioned nearly one 
synthesized beam away from the galaxy center and the 
total flux doubles on a timescale as short as 6 hours (see 
Fig. [8]). Typical scintillation effects at 325 MHz only 
create about 3% source variability on similar timescales. 
Therefore, we consider scintillation effects an unlikely 
source of the detected emission. 

To interpret the observation threshold for astrophys- 
ical sources implied by our transient detection, we ex- 
press measured radio flux density S in terms of a limiting 
source brightness temperature Tb ■ In the Rayleigh- Jeans 
limit, this can be written as 



T B = 



2k B v 2 VL s ' 



(4) 



where v is the observing frequency and Q s is the source 
solid angle. Expressing the source solid angle in terms of 
a characteristic size r and distance D, and inserting val- 
ues relevant to this search, Equation 2] can be rewritten 



T B - 1 x 10 8 K 



S 
mjy 



R, 



(5) 



where i? Q is the solar radius. The range of astrophysical 
sources probed by our transient detection is illustrated 
in Figure [TU1 

Additional limitations to the transient progenitor can 
be inferred by considering the temporal characteristics 
illustrated in Figure [H We observe a ~ 1 mjy increase 
in radio flux in 6 hours, implying an upper limit to the 
source size of c x 6 hr ~ 42 AU. Using Equation \E\ in- 
coherent emitters (approximate temperature boundary 
1 x 10 12 K) of this size scale are detectable out to 900 kpc. 
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We observe no evidence of variability on timescales 
less than one hour. This excludes sources with short 
tem poral scales, such as pulsars and related ph enom- 
ena ()Staelin fc ReifensteinH 19681 : iBhat et al.ll2007T) CMI 
burst emission from ultra cool dwarfs (B erger et al.ll2001fc 
IQsten et al.l [200 9), and exoplanets (Lazio ct al. 2004). 
Transient J103916. 2+585124 also appears to decay on 
a time scale less than one day. This may be ex- 
plained by the marginal detection of radio emission 
from supernovae. S pectral luminosity measurements by 
iWeiler eta l. (2002) indicate that radio supernovae may 
be detected by this transient search at distances exceed- 
ing 20 Mpc. 

We further consider the likelihood that our transient 
detection is due to coherent emission from a stellar flare . 
Observations of flare star EV Lac by Osten et al. (2005) 
revealed intense emission exceeding 10 13 K with temporal 
structure similar to that observed for J103916. 2+585124, 
i.e. minute-scale rise time and a decay time of hours. 
Further, flare stars are common in the solar vicinity, 
with a catalog of nea rby UV Cet-type flare stars by 
iGershberg et al.l ([1999) indicating 463 objects within 50 
pc. Outbursts exceeding 2 x 10 14 K with physical size 
scales equal to one solar radii should be detected in our 
observations for sources within 1 kpc. While there is 
no published stellar counterpart, our astronomical de- 
tection distance is roughly twice the reported 2MASS 
non-detection limit for EV Lac type stars. 

5. SUMMARY 

We present the results from a 325 MHz variability and 
transient search using data from the VLA archive. The 
search revealed multiple variable sources well described 
by interstellar scintillation of extragalactic sources and 
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